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Abstract 

Mode-locking regions (resonance tongues) formed by border-collision bifurcations of 
piecewise-smooth, continuous maps commonly exhibit a distinctive sausage-like geom- 
etry with pinch points called "shrinking points" . In this paper we extend our unfolding 
of the piecewise-linear case [Nonlinearity, 22(5):1123-1144, 2009] to show how shrinking 
points are destroyed by nonlinearity. We obtain a codimension-three unfolding of this 
shrinking point bifurcation for iV-dimensional maps. We show that the destruction 
of the shrinking points generically occurs by the creation of a curve of saddle-node 
bifurcations that smooth one boundary of the sausage, leaving a kink in the other 
boundary. 



1 Introduction 

Piecewise-smooth systems are used to model a vast range of physical systems involving 
nonsmooth behavior pfl [2j [3j H] . In this paper we study piecewise-smooth, continuous maps, 
i.e. 

x i+ i = F(xi) , (1) 

*DJWS acknowledges support from an NSERC Discovery Grant. JDM acknowledges support from NSF 
grant DMS-0707659. 
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where G M, N and F is everywhere continuous but nondifferentiable on codimension-one 
surfaces in R^ called switching manifolds. Such maps arise as Poincare maps of Filippov 
systems near grazing-sliding and corner collisions E] or of some time-dependent piecewise- 
smooth flows [7J. They are often used as mathematical models of various discrete, nonsmooth 
systems, see e.g. [8]. 

A fundamental and unique bifurcation of piecewise-smooth, continuous maps results from 
the collision of a fixed point with a switching manifold; it is known as a border- collision 
bifurcation. Except in degenerate border-collision bifurcation may be classified as 

either a border- collision fold at which two fixed points collide and annihilate, or a border- 
collision persistence at which a single fixed point "crosses" the switching manifold [TJ |9] . Note 
that the collision of one point of a periodic solution of ([Q) with a switching manifold is also 
a border- collision bifurcation for the n th iterate of ([T]) [HEED]. Though complicated dynamics 
may be born in border-collision bifurcations, in this paper we study only the creation of 
periodic solutions. 

We assume that the derivatives of the smooth components of (pQ), D X F, are locally 
bounded (we exclude from consideration, for instance, square-root type maps [HEEE!). Then, 
generic border-collision bifurcations of (JTJ may be described by piecewise-linear maps. More 
precisely, structurally-stable dynamics of a piecewise-smooth, continuous map near a border- 
collision bifurcation are described by the piecewise-linear, series expansion about the bifur- 
cation. A consequence is that, to lowest order, the structurally-stable invariant sets created 
at border-collision bifurcations grow linearly as the bifurcation parameter varies. 

One example is the two-dimensional, piecewise-smooth, continuous map 
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where 



s = eJx (the first component of the vector x & R 2 ) , 
< 7"l,sr < 1, < u\_, wr < ~, fi e R is assumed to be small, and the functions g L (x) 
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and 

x) contain the terms of f ^ that are nonlinear in x. The map ([2]) is continuous only if 

whenever s = 0, i.e., on the switching manifold. 
A border-collision bifurcation for ([2]) occurs at the origin when /x = 0. If this bifurcation 
is nondegenerate, the local dynamics are independent of the nonlinear components, g L and 
9 R - 

The piecewise-linear version of (jSJ), i.e., with 



g L (x) = g R (x) 








(3) 



is the canonical, piecewise-linear form for a border-collision analogue to a smooth Neimark- 
Sacker bifurcation. Indeed the fixed point for fi < has a pair of complex multipliers 
(A± = r L e ±2?rla;L ) inside the unit circle that "jump" at /j = outside the unit circle (A± = 
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— e ±27rl£jR ) for ji > 0. Depending upon the precise choice of parameter values, this border- 
collision bifurcation may generate periodic, quasiperiodic, or chaotic solutions as well as 
combinations of these [T21 EH [EU EE] • This class of border-collision bifurcations has been 
seen in models of DC/DC power converters [7] and optimization in economics [H IT6] . 
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Figure 1: Panel A shows a bifurcation diagram of ([2]) with ([3]) when r\_ = 0.2, sr = 0.95 
and Co>i_ = wr = 0.287. When /i < the unique fixed point is attracting (the solid line) and 
when fi > it is repelling (the dotted line) and a pair of period-seven orbits exist. The 
solid lines for /i > in panel A show the attracting 7-cycle (two pairs of points have similar 
s- values). Panel B shows a phase portrait when \i — 1. The dotted vertical line is the 
switching manifold. The attracting [saddle] 7-cycle is indicated by triangles [circles]. The 
unstable manifold of the saddle forms an invariant circle (red); the curves forming its stable 
manifold (blue) intersect at the repelling fixed point (square). 



Fig. [I]- A shows an example of a bifurcation diagram for ([2]) with when a pair of period- 
seven orbits (7-cycles) is created at fi = 0. These two orbits are shown in panel B; one is 
attracting and the other is a saddle. The unstable manifold of the saddle forms an invariant 
circle that contains the attracting orbit. This border-collision bifurcation is nondegenerate; 
indeed, if we were to relax ([3]) then there is an e > such that whenever < /i < e the 
dynamics are topologically equivalent to Fig. [T]-B near the origin. An example is shown in 
Fig. [2] where 

(4) 

For this nonlinear system, there is also a saddle-node pair of period-seven orbits created at 
the border-collision bifurcation. The stable 7-cycle is attracting up to fi ~ 1.202. 

Note that piecewise-linear maps are particularly straightforward to analyze because any 
periodic orbit is the solution to a linear system. Furthermore, linearity implies that if X is an 
invariant set of f ^ then AX is an invariant set of for any A > 0. Consequently it suffices 
to consider /i = —1, 0, 1. 

A two-parameter bifurcation diagram of the piecewise-linear case of ()2|) is shown in Fig. [3] 
for fj, = 1. The colored regions are resonance (or Arnold) tongues within which there is an 
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Figure 2: A bifurcation diagram showing attracting orbits of (j2J) with nonlinearity @ at the 
same parameter values in Fig. [TJ When \i is small the dynamical behavior near the origin 
is topologically equivalent to the linear case. However, at [A ~ 1.202, the attracting 7-cycle 
undergoes border-collision beyond which there exists a complicated attracting set. 



attracting periodic solution. Since the piecewise-linear case of the two-dimensional map (j2J) 
has a unique fixed point, we may define a rotation number for orbits as the average change 
in angle per iteration about the fixed point (T2J [171 US EE]- Consequently the resonance 
tongues can be labeled by the rotation number, m/n, of the corresponding periodic solution. 
The orbits shown in Fig. CD lie in the 2/7-tongue; this tongue intersects the Sr — 1 line at 
ojr = 2/7 « 0.2857. 

The majority of the resonance tongues in Fig. [3] exhibit a structure that is often likened 
to a string of sausages. This structure was first observed in a one-dimensional sawtooth map 
|21j . and has since been described in higher dimensional maps such as (|2J), see for example 
[HI 1221 El DSL As in |2"Tj . we refer to points where resonance tongues have zero width as 
shrinking pomto 

An unfolding of shrinking points in piecewise-linear, continuous maps of arbitrary dimen- 
sion was performed in [10J. There it was shown, upon imposing reasonable nondegeneracy 
assumptions, that any two-dimensional slice of parameter space in the neighborhood of 
a shrinking point will resemble Fig. HJ shrinking points are codimension-two phenomena 
of piecewise-linear, continuous maps. In particular, near a shrinking point, the resonance 
tongue is locally a two-dimensional cone bounded by four curves. These boundary curves are 
pairwise tangent at the shrinking point so the cone boundaries are locally C 1 . In the interior 
of the cone a primary n-cycle exists; it has some number of points, say I, located, "left", of 
the switching manifold (1 = 2 for Fig. H]). This orbit collides and annihilates with another 
n-cycle in a border-collision fold bifurcation on the four boundary curves of the cone. This 
secondary periodic solution has I — 1 points to the left of the switching manifold within one 
half of the cone and / + 1 points to the left of the switching manifold within the other half 

^dn this paper we will only consider the case of "non-terminating" shrinking points defined in 110] . Thus 
we do not consider the ends of resonance tongues such as those at Sr = 1 in Fig. [3] 
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Figure 3: Resonance tongues of fl2]) with ([3]) for ty = 0.2 and \i — \. Each colored region 
corresponds to the existence of a periodic solution with period indicated in the color bar. The 
rotation number in each resonance tongue is equal to the value of cjr at which it emanates 
from the line, sr = 1. White regions correspond to no observed attractor; black regions 
correspond to the existence of a bounded forward orbit that is either aperiodic or has period 
larger than 30. Similar figures at different parameter values are given in [121 120] . 



of the cone. On each of the four boundary curves a different point of the primary orbit lies 
on the switching manifold. 

These results generalize to the nonlinear case in the follwing sense: regardless of the 
nonlinearites g L and g R , the two-dimensional bifurcation structure shown in Fig. [3] will be 
essentially unchanged when /x is small enough and it will limit to the linear picture as 
/i — » + . However, when nonlinear terms are present the sausage structure is not preserved 
as /i increases; indeed, the shrinking points break apart as shown in Fig. |5j Consequently 
when nonlinear terms are present the phenomenon is codimension-three; we refer this scenario 
as a generalized shrinking point. 

Fig. [5] suggests that the break-up of shrinking points occurs in a regular fashion. The 
resonance tongues develop a nonzero width that increases with [i. For \i > the right 
boundary of each resonance tongue appears to be smooth whereas each left boundary retains 
the kink that appeared at the shrinking point. It is interesting that, in contrast to the case 
for smooth maps, we observe no "strong" resonance behavior when n < 5. 

The purpose of this paper is to determine the generic unfolding of generalized shrinking 
points for a piecewise-smooth, continuous map of arbitrary dimension. Here we summarize 
our results. Suppose that a border-collision bifurcation occurs when a parameter /x is zero, 
that resonance arises for \x > 0, and that in the limit \x — > + , a two-parameter bifurcation 
diagram exhibits a shrinking point with its four boundary curves. For small /i > the 
four boundary curves maintain a common intersection point, O, but each curve is no longer 
necessarily tangent to the opposite curve at this point. 

We will show that, under reasonable nondegeneracy assumptions, there exists a new 
bifurcation locus for each small enough fixed \i > 0. This locus is a curve of saddle- node 
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Figure 4: A magnification of the 2/7 resonance tongue in Fig. [3] near a shrinking point. 
Orbits that exist in the resonance tongue are shown in inserted phase portraits (not to scale) 
as the open circles and triangles. Above the shrinking point, these have / = 1 and I = 2 
(I represents the number of points that lie left of the switching manifold) and below the 
shrinking point, / = 2 and / = 3. On the four boundary curves the two coexisting periodic 
solutions collide and annihilate in border-collision fold bifurcations. On each of the four 
boundary curves a different point of the primary orbit (with 1 = 2), lies on the switching 
manifold. These orbits are sketched in phase portraits with filled circles. Compare this to 
panel B of Fig. CEJ an actual phase portrait at (ur, sr) = (0.287, 0.95). 



bifurcations of the primary n-cycle that is tangent to one boundary curve at a point A, and 
to an adjacent boundary curve on the other side of the shrinking point at B, see Fig. |6j This 
locus is smooth and collapses to the shrinking point as — > + . Two n-cycles that have 
the same itinerary as the primary n-cycle exist in a triangular region bounded by O, A and 
B. If we let Q\ denote the angle made at O between the two border-collision curves across 
the region AOB, see Fig. [6j and 82 is the opposing angle made at O between the other two 
border-collision curves, then for small u > 0, 61 < 8 2 . 

A formal statement that includes these results is given in Thm. [10] in £j6j 
Here we summarize the remainder of the paper. In the following section we present the N- 
dimensional map (EJ) that describes an arbitrary border-collision bifurcation. Concepts from 
symbolic dynamics that are invaluable for describing periodic solutions of are introduced 
in Key formulas for periodic solutions are obtained in §HJ Subsequently we impose the 
assumption that <Q is piecewise-C^ ; this allows us to derive useful series expansions. Section 
[5] is devoted to defining generalized shrinking points as a codimension-three scenario. Finally 
in §0 we unfold these points leading to Thm. \IU\ The proofs of the theorem and Lemma [9] 
are given in an appendix. 

Throughout the paper we use 0{k) [o{k)\ to denote terms that are order k or larger 
[larger than order k] in all variables and parameters of a given expression. 



6 




Figure 5: Resonance tongues of ^2} with (jlj) when r\_ = 0.2 for three different values of 
/i. The uppermost plot is a magnification of Fig. |3j The lower two plots were computed 
numerically by estimating the eventual period of the forward orbit of the origin. 

2 Generic border-collision bifurcations 

We restrict our attention to local dynamics of border-collision bifurcations, and thus study 
a piecewise-smooth series expansion of (JT|) about an arbitrary border-collision bifurcation. 
Throughout this paper we will assume that there is a single, smooth switching manifold in 
a neighborhood of the border-collision bifurcation. This is typically the case in models since 
switching manifolds are usually defined by some simple physical constraint. For simplicity we 
assume that a coordinate transformation has been made so that the switching manifold cor- 
responds to the vanishing of the first component of x (see in particular [23] for details of 
such a transformation) and — to avoid the use of subscripts for components — we introduce 
the new variable 

s = ejx . (5) 
A general piecewise-smooth map then takes the form 

*m = /(*,;£) = { f f R^\ ' (6) 

where £ is a vector of parameters. The switching manifold partitions phase space into two 
regions: the left half-space (where s < 0) and the right half-space (where s > 0). 
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Figure 6: Bifurcation sets of ([2]) with (j3J) for tl = 0.2 and three different /i values. Panel A 
is identical to Fig. HI In panels B and C the locus of classical saddle-node bifurcations of the 
primary orbit is indicated by a red curve. Solid [dotted] curves correspond to border-collision 
fold bifurcations [border-collision persistence]. 



We assume a border-collision bifurcation of a fixed point occurs at the origin when a 
parameter /i is zero so that the functions f L and / R have series expansions 

f J (x t ; = ^6(0 + Aj(Z) Xi + g J (x i; , (7) 

where J = L, R, \i is the first component the parameters £, Aj is an iV x iV matrix, b e Mr. 
The functions g J contain only terms that are nonlinear in x; that is, g J (xi] £) = o(|xj|) (this, 

3 

for example, could include terms of order \Xi\2 that arise naturally in Poincare maps relating 
to sliding bifurcations). By continuity of fl6]), A\_ and Ar are identical in their last N — 1 
columns and g L = g R whenever s = 0. 

If Aj(£) does not have an eigenvalue 1, the half- map f J has a unique fixed point near 
the origin given explicitly by 

x* J (£) = (I - MtVXt) A* + o(Ai). (8) 

As seems to have been first noted by Feigin, see [9], a convenient expression for the first 
component of the vector, x* J , is obtained with adjugate matrices [2^, [25]: 



adj(X)X = det(X)/ , for any N x N matrix X. (9) 

The point is that since A\_ and Ar are identical in their last N — 1 columns, adj(J — A\_) and 
adj(J — Ar) share the same first row: 

g T (0 = e7adj(J - A L (£)) = e]"adj(J - A R (0) • (10) 

Consequently, multiplication of (IH]) by e]" on the left implies that the first component of the 
fixed point x* J satisfies the useful formula 



det(I-Aj(0) 
8 



// + o(//). (11) 

/i=0 



In particular we learn from (11 ip that both fixed points, if they exist, move away from the 
switching manifold linearly with respect to /x if and only if £> T (£)&(£) 1^=0 7^ 0. This condition 
is a nondegeneracy condition for the border-collision bifurcation. Under this assumption, the 
bifurcation is a border-collision fold bifurcation if det(I — A|_(0)) and det(J — Ar(0)) have 
opposite signs, and border-collision persistence if they have the same sign j9]. 

3 Symbolic dynamics 

It is common in the study of piecewise-smooth systems for symbolic methods to be highly 
beneficial. In this paper we consider bi-infinite sequences, S, constructed from the binary 
alphabet {L, R}. In order to unfold shrinking points, we find it necessary to consider only 
what we have termed rotational symbol sequences. In (TO] we defined these as particular 
finite collections. Instead of repeating this definition we provide here a definition that is 
more versatile in that it extends naturally to nonperiodic sequences and has been described 
elsewhere. Furthermore, to be consistent with combinatorics literature sequences are always 
assumed to contain infinitely many elements (unlike in [TO]). 

Definition 1. For a,j3 G [0,1), let S[a,/3], be the bi-infinite symbol sequence with i th 
element 

S[a,P)i = 



ia mod 1 G [0, (3) 
iot mod 1 G [(3, 1) 



for all i G Z. 




Figure 7: A geometric portrayal of Def. [TJ First cut a circle with a vertical line that subtends 
an angle 2tc(3 as shown. Each real number, 0, is then represented by a point 2ir(j) radians 
clockwise from the lower intersection of the circle with the vertical line. Then Si[a, 0\ = L 
whenever = ia is located to the left of the vertical line, and Si[a, 0\ = R otherwise. In 
addition <So[a, 0\ is always L. 
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A visual representation of these sequences is provided by Fig. [0 These sequences seem 
to have first been studied by Slater [2SJ [27| • The sequences for which a = (3 £ Q have been 
well studied and are known as rotation sequences; they are equivalent to Sturmian sequences 
[28] . which are traditionally defined from a combinatorics viewpoint [29]. Some discussion of 
the case a = (3 G Q is given in [30]. A related approach is to consider arithmetic sequences 
on Z/nZ for some n G N [31] . Similar sequences for rotational orbits of the Henon map are 
described in [32J. 

Definition 2 (Rotational Symbol Sequence). For each l,m,n G N with l,m < n and 
gcd(m, n) = 1, <S{^, -] is a rotational symbol sequence. 

Any periodic sequence, such as S[^, j-}, is generated by repeated copies of a finite col- 
lection from {L, R}, termed a word. For instance for the word W = LLR, the generated 
sequence is S = . . . LLRLLRLLR. . .. Here we let W[l,m, n] denote the word comprised of 
the i — 0, 1, . . . , n — 1 elements of S\— , -1. 

As an example, let us compute <S[I, |]. Here a = |, thus the numbers ia mod 1 of Def. [1] 
for i = 0, 1, . . . , n - 1 are 0, f , f , f , f , f , f , and /3 = f , therefore W[3, 2, 7] = LLRRLRR. 

Throughout this paper will we use the symbol d to denote the multiplicative inverse of 
m modulo n, for example d = 4 when m = 2 and n = 7 as above. Then 



where id is taken modulo n. For clarity, throughout this paper we omit "modn" where it is 
clear modulo arithmetic is being used. 

Given a word W, we let denote the i th left cyclic permutation of W and W l denote 
the word that differs from W in only the i th element. For example ifW = LRLRR then 

= LRLLR and W (2) = LRRLR. If S is the sequence generated by W we let S w and 
denote the sequences generated by and W\ respectively. 

A key property of symbol sequences that is verified in [10] is 



4 Periodic solutions 

For any symbol sequence S, the iterates of a point xo G under the two half-maps of ([6]) 
in the order determined by S are 



In general this may be different from iterating xq under the map (J6]). However, if the sequence 
{xi} satisfies the admissibility condition: 




(12) 




(13) 



(14) 




L, whenever Sj < 
R, whenever Sj > 



(15) 
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for every i, then {x{\ coincides with the orbit of Xq under ([6]). Notice if = there is no 
restriction on «Sj. When (TT5]) holds for every i, {xi} is admissible, otherwise it is virtual. 

When S has period n, the periodic orbits with this sequence are admissible fixed points 
of the map 

f = f- 1 o-..of°. (16) 
Some straightforward algebra leads to 



f{x; = P 5 (0&(0 » + o{p) + (m s (0 + o(fi )) x + o(x 



(17) 



where 

M s = As^.-.Aso, (18) 
Ps = I + A Sn _ 1 +A Sn _ 1 A Sn _ 2 + --- + A Sn _ 1 ...A Sl . (19) 

Notice Ps is independent of S , thus 

Ps = P s o • (20) 

We now present six lemmas relating to periodic solutions that are useful for analyzing 
shrinking points. A consequence of the following lemma is that M s and M s0 - are identical 
in their last N — 1 columns. 



Lemma 1. For any N x N matrix, X , XAr = XAi + (ej, for some ( G 



Proof. Since Ai and Ar are identical in all but possibly their first columns, we may write 
^4r = A\_ + C e i" f° r some £ G M. N . This proves the result with ( = X(. □ 

If Xq is a fixed point of f s then the n points {a;o, . . . , £ n _i} (where Xi is defined by f|T4l) ) 
describe a periodic solution (which may be virtual) that we will refer to as an S-cycle. The 
following two lemmas are generalizations of those given in |10j . 



Lemma 2. Suppose {x^ is an S-cycle and Xj lies on the switching manifold, for some j. 
Then {x{\ is also an S ] -cycle. 

Proof. By continuity: f L (xj]l;) = f R (xj]l;), hence there is no restriction on the j th element 
of S. □ 

Lemma 3. Suppose {x{\ is an S-cycle. Then for any j , det (/ — D x f s<J \xj; £)) = det (/ — 
D x f(x ;0)- ' 

Proof. By the chain rule the Jacobian D x f s (xo~, £) may be written as a product of matrices: 

o 

D x f(x ;0= II I >-.f (■<:■. 0- 



i=n—l 



The spectrum of this product is unchanged if the iV matrices are multiplied in an order that 
differs only cyclically from this one (refer to the proof of Lemma 4 of [10] ) , which proves the 
result. □ 
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Whenever Ms(£) does not have an eigenvalue 1, the implicit function theorem implies 
that f s has a unique fixed point near the origin, call it x* s . Using (1171) we obtain 

x* s (0 = (I - MsiOr'PsiOKO V + o(ji) • (21) 

/j,=0 

We may derive a formula for the first component of x* s , denoted s* s , in the same spirit as 
(TTTj) for fixed points of (jSJ). Let 

^(0 = e[adj(/-M 5 (0), (22) 

then 



det(/-M s (0) 



M + (23) 

At=0 

As in [20] we use two lemmas to derive a convenient formula for s* s . 



Lemma 4. The matrices Ps{I — As ) an d I ~ Ms can differ in only their first columns. 
Proof. Using (jig} and flI5}. 



Ps(/-As ) = {I + A Sn _ 1 +As n _ 1 As n _ 2 + --- + As n _ 1 ...A Sl ){I-As ) 

expand and group terms differently: 
= I-M s + (As^ - A So ) + As n _M Sn _ 2 -A So ) + --- 

+ A Sn _ 1 . . . A S2 (A Sl - A So ] 

apply Lemma [H 
= I-M s + (n-iej + C- 2 e]" + • • • + CieJ 
where each Q e IR N , 



'n-1 



I-M s + [Y,Ci Ui ■ 



,i=i 



□ 



Lemma 5. f?J-Ps = det(Ps)o T ; where Qg is given by ((Hj) and o T is given by [TU) . 
Proof. By Lemma H] we have 

e]adj(P 5 (/-4s )) = ejad](l - M s ) = 

e]"adj(J - A So )&d](P s ) = Qs (since adj(XF) = adj(F)adj(X) for any X, Y) 

£? T adj(Ps) = g T s by<m 

det(P s )g T = g T s P s by ® 



□ 
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If I — Ms(£) is nonsingular, by Lemma [5] and (1231) . 



s * S(e) = det(P g (Q) T 



fi + o(/i) . (24) 

At=0 



This expression relates the linear component of s* 5 (£) simply in terms of g T b (which appears 
in the fixed point equation (|lip ) and the determinants of P$ and I — Ms- (Feigin's result 
concerning the creation of 2-cycles at border-collision bifurcations (see [H]) follows from (|24|) 
by substituting S = LR and S = RL). 

We finish this section with an important lemma that is most simply stated for ([6]) in the 
absence of nonlinear terms, for then all o(/i) terms given above vanish. Though this result 
is proved in (JO] , with the use of Lemma [5] we are now able to provide a pithier proof. 

Lemma 6. Suppose the map is piecewise-linear, that is g L = g R = 0, and assume /i ^ 
and g T b ^ 0. 

i) If I — Ms is nonsingular, then the unique fixed point of f s , x* s , given by WT\) . lies on 
the switching manifold if and only if Ps is singular. 

ii) If Ps is nonsingular, then f s has a fixed point if and only if I — Ms is nonsingular. 



Proof. Part (JI]) follows immediately from fl24|) . If / — Ms is nonsingular, part (JD} is trivial 
by fl2T|) . Suppose / — Ms is singular and suppose for a contradiction f s has a fixed point 
x* s . Then by f[T7|) we have 

(/ - M s )x* s = P s bfi . 
Multiplication of this by g~s fl22l) on the left yields 

det(J - M s )s* s = det(P s )g T bn , 

where we have also used (JH]) and Lemma This provides a contradiction because the left 
hand side of the previous equation is zero, whereas by assumption the right hand side is 
nonzero. □ 



5 Generalized shrinking points 

At this stage we find it useful to impose the extra assumption that the map under 
investigation here is piecewise-C^ , for some K G N. This allows us to derive series expansions 
of smooth components of the map and iterates of the map. In particular this assumption 
allows us to apply the center manifold theorem necessary for proving the existence of saddle- 
node bifurcations in §|6j 

In order to unfold a generalized shrinking point, we must first give a precise definition 
of such a point. We use the results of the previous section to write down assumptions that 
guarantee the existence of a periodic solution with two points on the switching manifold. 
As in [10], it is useful to assume that this periodic solution is admissible. To state this 
assumption we need to renormalize the map (J6]). 
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Recall that the map (J7|) depends upon an arbitrary vector of parameters £, and that \i 
denotes the first component of this vector. Scaling x by \x gives a new map h defined through 

f J (jiz\g) =nh J (z\g) , 

where z G M. N . Note that when f J is C K then h J is C K ~ 1 , and using (j7j), it has the expansion 

h J (z;t)=b{t)+Mt)* + 0(lM). 
For /i > 0, the renormalized map for z is then 

z n+1 - h(z n , S)-< ^ , (25) 



where 



u = ejz 



Whenever /i > 0, if x — /j,z then = f/,h(z',£). For any 5-cycle (12TT) . we can also 

let = so that z* s (£) is C K_1 and is a fixed point of h s = /i 5 "- 1 o • ■ • o h s °. 

Since (1251) is a "blow-up" of phase space, points z G M. N are not necessarily near the origin 
when fi is small. The renormalization effectively transfers the /^-dependence of the piecewise- 
smooth map from the constant term to the nonlinear terms and this scaling is often helpful 
in the analysis. Note that (1251) is has nontrivial dynamics for /x = 0; indeed in this case it is 
piecewise-linear and identical to the map studied in [10J. 

We now use Lemma M to guarantee the existence of a periodic solution with two points 
on the switching manifold in terms of singularity of matrices, P s w (jT9l) . 

Definition 3 (Generalized Shrinking Point). Consider the map © with N >2 and suppose 
f? T (0)6(0) 7^ 0. Let S = S[— , ~] be a rotational symbol sequence with 2 < I < n — 2. Suppose 

-Ps(O) and Psw-iw (0) are singular (the singularity condition). 

Let 

S = S 5 , (26) 
S = S m , (27) 

and assume I — Mg(0) and / — M§(0) are nonsingular. Let {£«(£)} be the unique iS-cycle 
near the origin (xo(0 is given by (12TT)). Let 

Vi = ^(0) , (28) 

where = ^ the or bit {yi} is admissible then we say that (jHJ) is at a generalized 

shrinking point when £ = 0. 

The sequences S ( 12 6p and 5 (T27|) are rotational with one less and one more L than 
S = S[22, 1], respectively (specifically, £ = , and <S = <S[f , *±*]). The periodic 

solution is fundamental to the shrinking point. As stated in the following lemma, it has 
two points on the switching manifold. Moreover if {£«(£)} denotes the unique iS-cycle near 
the origin, then ij(0) = Zi(0) = yi. We let 

U = ejyi . (29) 
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Lemma 7. Suppose (0|) is at a generalized shrinking point when £ = 0. Then, 

i) t = t ld = 0; 

H) td,t(i-i)d < 0, t~d,t(i+i)d > 0; 
Hi) {yi\ has minimal period n; 

See [10] for a proof. Lemma [7J essentially states that the orbit {yi} appears as in Fig. [HJ 
A consequence of Lemma [7J is that several important matrices are singular. To see this, 
first note that the point y is a fixed point of h s (y;0). By Lemma [2] and Lemma [TJI]) , y 
is also a fixed point of h s ° ld (y;0). Using fTl2|) . S old = S^ d \ and so yd is a fixed point of 
h s (y; 0). The points y and y d are distinct (by Lemma fTT pTi]) and since is admissible), in 
other words there are multiple 5-cycles. Hence the matrix / — M$(0) must be singular and 
consequently each P S (t) (0) is singular by Lemma EJJii]) producing the following result (given 
also in [TO]): 

Corollary 8. Suppose (OJ) is at a generalized shrinking point when £ = 0. Then, 
i) I — Ms(0) is singular; 

ii) P S (j)(0) is singular, for all i. 

For reader convenience let us briefly summarize symbols used: 

x G M. N , s = ejx, 

yi = Zj(0) = ^(0) , t = e[y , 

fiz = x , it = e x ^ . 



'a-i)d 



'd+l)d 



'-d 



Figure 8: The orbit {y^} as described by Lemma [7J The points y$ and yid lie on the switching 
manifold, s = 0. 
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6 Unfolding generalized shrinking points 



We begin by performing a change of coordinates, similar to that in [lOj, such that, locally, 
two boundaries of the associated resonance tongue lie on coordinate planes. We are given 
that tto(0 and uid(Q are C K ~ 1 and Uo(0) = w/d(0) = (Lemma W($))- Since a generalized 
shrinking point is a codimension-three phenomenon, we assume there are three bifurcation 
parameters 



As long as the matrix 



diiQ dii ' 



drj 
duid 



dv 
duid 



drj dv . 



is nonsingular, we may perform a nonlinear coordinate change such that 

fco(0 = 77(1 + 0(1)), 
MO = v(l + 0(l)). 



(30) 
(31) 



Consequently, on the coordinate plane rj = 0, the point Xq, of the 5-cycle lies on the switching 
manifold. According to Lemma [21 the 5-cycle is also an S° = 5-cycle here. Similarly 
on v = 0, x id lies on the switching manifold. By Lemma [2J here the 5-cycle is also an 
gid _ S^~ d >-cyc\e (this equality follows simply from ( I12p ). Along the /i-axis (in three- 
dimensional parameter space) both xq and xid lie on the switching manifold so here the 
iS-cycle is also an iS-cycle. Hence on the /i-axis the orbits {x^)} and {x^)} are identical. 

In order for £ to properly unfold a generalized shrinking point we need a nondegeneracy 
condition on the nonlinear terms of ([6]) that guarantees the shrinking point breaks apart in 
the usual manner as /i increases from 0. Along the /i-axis, xq is a fixed point of f s and when 
fi = it has an associated multiplier of 1. We will assume that the algebraic multiplicity of 
this multiplier is one. An appropriate condition on the nonlinear terms is that this multiplier 
varies linearly (to lowest order) with respect to /j, for small /j. That is, we require 



|-det (I-D x f(x (0;0) 



7^0 



(32) 



We are free to scale the parameter /i before performing the analysis. However, we find it is 
most instructive to merely fix the sign of fx in a way that ensures resonance arises for fi > 0. 
It turns out that the following assumption achieves this effect: 



3gn (J- det (/ - D x f(x (0; 0) \) = sgn(det(J 



(33) 



Our unfolding theorem details the existence of admissible periodic solutions in regions of 
three-dimensional parameter space that are bounded by six different surfaces. As a conse- 
quence of the choice (|3ll and (I3"T|) . three of these surfaces are simply the coordinate planes. 
The following lemma gives series expansions of the remaining three surfaces. 
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Lemma 9. Suppose the piecewise-C K map (TJJ) is at a generalized shrinking point when £ = 
and that K > 4. Assume that the only eigenvalue of Mg(0) on the unit circle is 1 and that 
it has algebraic multiplicity one, and that (fffiD)) . I[JJ\) and (f^j /io/d Let 



5 = det(J - Mg(0)) , 
5 = det(J - M^(0)) , 

ft 

h : 



— det(I-D x f(x (0;0) 



(34) 
(35) 

(36) 



(which are all nonzero due to assumptions in the definition of a generalized shrinking point). 
Then, 



i) 



S 



tdt(i-i)d 



8 t-dt(l+l)d 

H) uid(£) = on a C K ~ 1 surface, r\ = v) — v 



kot d 



t d 



5t(i +1 ) d t{i-i)dt(i+i)d 



u + 0(2) ; 



Hi) uq(0 = on a C K 1 surface, v = 2 (/i,^) = r\ 



iv) there is a C function given by 



k t 



(l-l)d 



t 



st- d 



ix- 



(Z-l)d 
tdt-d 



s 



td 



t 



(Z-l)d 



A + ( 2 ) 
0*d 



(37) 



suc/i £/ia£ /or /i > classical saddle-node bifurcations of S -cycles occur when A(£) = 
(though are not necessarily admissible); 

v) if n > 0, then A(£) < only if rj < and z/ > ^{fi, v) > moreover A(£) = a/ong 
(//, 0, Ci(/i)) and (//, C 2 (/i), 02 (A*, C2G"))) /or C^- 1 functions (i and Q 2 , 



CM 

CM 



kot(i-i)d 



5 



li + 0{ii 2 ) , 



&0*d 



fi + 0{n 2 ) . 



(38) 
(39) 



See Appendix |X] for a proof. It is useful to consider the nonsingular, linear coordinate 
change: 

~ ko 1 1 
5 



/' 

v 








td t(l-l)d 
1 1 



t d 
1 

id 



t(i-i)d 



t(i-i)d 



V 



(40) 
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because substitution of (j4"0]) into (|37j) produces 

A(£) = /i 3 + tp - P 3 + o(2) . (41) 

In this alternative coordinate system, saddle-node bifurcations occur approximately on a 
cone. All six surfaces are sketched for fi > in Fig. [9j The surface A(£) = intersects the 
plane r\ = tangentially along the curve (/^, 0, Ci(aO) and the surface v = 2 v) tangentially 
along the curve (//, ( 2 (a0, 02 (a*, C2 ()"))) • 

Different two-dimensional slices of parameter space through Fig. [9] will produce vastly 
different bifurcation sets. Since slices defined by fixing the value of fi have been shown earlier 
(see Fig. E]), in Fig. Owe draw a plane at fixed \i > and show its curves of intersection 
with the nearby surfaces. This cross-section is shown again in Fig. [TOJ A close inspection 
of the formulas for (pi and 02 given Lemma [9] reveal a specific geometrical arrangement near 
the generalized shrinking point as follows. Since t d < 0, > (Lemma and 

sgn(fc ) = sgn(<5) (1551) . the coefficient for the \iv term of 0i(/i, v) is positive. Consequently, 
for small fi > the angle 9 2 in Fig. [10] is greater than 90°. Similarly the coefficient for the \ir\ 
term of 02 if-,v) is negative and so 81 < 90°. Any smooth distortion of Fig. [TU1 will preserve 
the property Q\ < 9 2 . The curves 77 = 0i(/i, v) and i/ = 02(/i, 77) bend left and down because 
the v 2 term of 0i(/i, z/) and the 7/ 2 term of 02 (/•£, ??) are both negative. The saddle- node locus 
is approximately a parabola. 

Finally, we present the main result. See Appendix [A] for a proof. 

Theorem 10. Suppose (TJ|) is at a generalized shrinking point when £ = and that K > 4. 
Assume that the only eigenvalue of Ms(0) on the unit circle is 1 and that it has algebraic 
multiplicity one. Assume 7/0 and are the only points of {y{\ that lie on the switching 
manifold. Assume we have KSfy) , l[3l\) and / fffffj) . Let 

m l = {£ g R 3 I fi > 0,7] > 0,1/ > 0} , 

^2 = e R 3 I /i > 0,77 < 0i(/i,z/),z/ < M^v)} , 

^ 3 = {£ e JR 3 I /i > 0, ( 2 {v) < 77 < 0, 2 (/i, rf) <u < Ci(/i), A(£) > 0} , 

/or the functions described in LemmalQ (see Fig. WO)) . 

Then, in ^i, S and S-cycles are admissible and collide in border- collision fold bifurcations 
on v = and 77 = for v > Ci (/•*); in ^2, S and S-cycles are admissible and collide in 
border- collision fold bifurcations on 77 = 0i(/i, z/) an<i v = 02(/x, 77) /or 77 < ^{t^), in ^3, two 
S-cycles are admissible and collide in classical saddle-node bifurcations where A(£) = 0; the 
boundaries 77 = and v = 2 (/i, 77) correspond to border- collision persistence. 



7 Conclusions 

We have studied resonance arising from an arbitrary border-collision bifurcation of a piecewise- 
smooth, continuous map. When a periodic solution created in a border-collision bifurcation is 
rotational, in the sense described in §3j the corresponding resonance tongue typically has 
a sausage-like geometry, see Fig. [3J Shrinking points break apart as parameters are varied 
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Figure 9: The three-dimensional unfolding of a generalized shrinking point for the piecewise- 
smooth, continuous map (J6]). The surface /i = corresponds to the border-collision bifur- 
cation of a fixed point. On each of the four surfaces, 77 = 0, v = 0, 77 = v) and 
v = (f)2{^,f]) (which share a mutual intersection with the /z-axis), one point of an 5-cycle 
lies on the switching manifold. The surface A(£) = corresponds to classical saddle-node 
bifurcations of an iS-cycle. This surface is approximately a cone, as made evident by transfor- 
mation to the alternative (/},?/, z>)-coordinate system, (HOl) . and tangentially intersects r\ = 
and v = 2 (/x,77) along the curves (/z, 0, CiXaO) an d (/i, (2(H), fain, (2^))), respectively. The 
plane fi = const, illustrates the cross-section shown in Fig. [TUJ 



to move away from the border-collision bifurcation, Fig. We have proved Thm. [TU] which 
details the manner by which shrinking point destruction occurs. The results of the theorem 
are in complete agreement with numerical results, Figs. [5] and |6j Theorem [10] does not pro- 
vide an understanding of global properties of resonance tongues, for instance the observation 
that the majority of, or perhaps all of, the kinks in the resonance tongues of Fig. appear 
on the left sides of the tongues. 

If the map ([6]) has an invariant circle that intersects the switching manifold at two points 
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*1 




Figure 10: A two-parameter bifurcation diagram for the map ([6]) near a generalized shrinking 
point for fixed, small /i > 0. As stated in Thm. [TOj an admissible 5-cycle coexists with, (i) an 
admissible 5-cycle in (ii) an admissible 5-cycle in \l/ 2 ; (hi) a second admissible 5-cycle 
in ty 3 . Along the solid curves the two coexisting, admissible, periodic solutions undergo a 
border-collision fold bifurcation except on the curve connecting the points v = Ci(a*) an d t] = 
C2(a*) which corresponds to a classical saddle- node bifurcation of the two periodic solutions. 
The boundary between ^ and ^3 and the boundary between ty 2 and ^3 correspond to 
border-collision persistence. 



and the restriction of the map to the circle is a homeomorphism, then any periodic solution 
on the circle will have a corresponding symbol sequence that is rotational. Our results do 
not apply to periodic solutions born in border-collision bifurcations that are non-rotational. 
However such periodic solutions seem to be, in some sense, less common [12]. Note that 
we make no requirement that (jSJ) is a homeomorphism or has an invariant circle, only that 
corresponding periodic solutions are rotational. 

A limitation of Thm. [10] is that it includes the assumption that the map ([6]) is piecewise- 
C . Poincare maps relating to sliding phenomena are generically piecewise-smooth, contin- 
uous with a |, 2, |, . . . type power expansion [5j and hence not apply to Thm. [TD] It seems 
reasonable that in this case a similar result with different scaling laws could apply. A major 
hurdle in the analysis is that the center manifold theorem (applied in the proof of Thm. [TO]) 
is not immediately applicable to non-integer power expansions. 

In our definition of the codimension-three generalized shrinking point (Def. [3]), we require 
that both matrices / — M§ and I — M§ are nonsingular at £ = 0. We then prove Lemma [9] 
by computing parameter-dependent series expansions (see Appendix [A]) . However part (E]) 
of Lemma [9] is independent of the parameters and it seems there should be a more direct 
proof of this result. This problem remains for future work. Related problems that remain to 
be fully understood include the persistence of invariant topological circles created at border- 
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collision bifurcations and the simultaneous occurrence of a border-collision bifurcation with 
a classical Neimark-S acker bifurcation. 



A Proofs for Section [6 



Proof of Lemma® We divide the proof into four major steps. In the first step we use the 
formulas for wo(0 (1301) and wm(£) (|3T1) to derive formulas for u and uid, enabling us to 
compute the border-collision boundaries r\ = 0i(«, v) and v = 2 (",??)• The method used 
is an extension of Lemma 3 of [10] to account for the nonlinear terms in ([6]). In step 2 we 
continue this methodology to derive identities relating coefficients, verify (JI|) of the lemma 
and refine the formulas for <pi and 02- in the third step we compute the one-dimensional 
center manifold of f s to identify saddle-node bifurcations on the surface A = 0. Lastly, in 
step 4 we study at the intersections of A = with rj = and v = fai^,?]) to complete the 
proof. For convenience we write 

det(J - M s (0, 77, v)) = k x -n + k 2 v + 0(2) . (42) 

Step 1: Derive formulas for uq and Uid and define <f>\ and <p2- 
Periodic solutions of with associated symbol sequences that differ by a single symbol 
may be related algebraically (as shown below). However the sequences S and S differ in two 
symbols, so in general we cannot effectively compare the S and iS-cycles. But if one of the 
four important quantities, Mo, u^, Uo and ui d is zero, then one of the two n-cycles is also an 
5-cycle and the two n-cycles may indeed be effectively related. The four curves which make 
up the edges of a resonance tongue near a shrinking point correspond to where these four 
quantities are zero. We separate this step of the proof by where different assumptions on 
the parameters are made. 

a) Suppose r] = 0. Then by ( !30l . uo = 0. Thus the iS-cycle, {xi}, is also an iS-cycle. 
So each Xi is a fixed point of / 5W . In particular, xu is a fixed point of f s(ld) ^ which in the 
renormalized frame ( 12 5 p says that zid is a fixed point of h s{u) : 

zi d (/i,0,u) = h s{ld \z ld {(i,0,v);[i,0,v) ■ 

Expressing h s<ld) as a Taylor series centered at yid (|2"g|) and evaluated at Z\d yields 

Zid(fJ., 0, u ) = hSild) (Vidl 0, v) + D z h s(ld) {yid] u, 0, v) (i M (u, 0, v) - y id ) +0(3) , (43) 

where the next term in the expansion is 0(//)0(|i;^(//, 0, v) — yid\ 2 ) which is 0(3) because 
zid(0, 0, 0) = y ld . 

Now, zid is a fixed point of h s(ld) (for any sufficiently small £), therefore 

z ld (p, 0, u) = h s{U) {y ld ; u, 0, v) + D z h s(U) (y ld ; fi, 0, v) (4*0", , z/ ) ~ Vid) + 0(3) , (44) 

where the right-hand side is the Taylor series of h sild) centered at y ld . But = 
(refer to $3§, thus whenever u = ejz = 0, we have h s(ld) (z; £) = h s(ld) (z; ^) (by continuity of 
(jSJ))- Since ti d = ejyu = (Lemma dj!])), we have 

h sm (y ld ;n,0,v) = h S(ld \y ld ;ti,0,v). (45) 
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Furthermore, although in general the last N—l columns of D z h s and D z h s are not equal, 
since h s< ~ ld) = hp {ld) on the switching manifold, the last N — l columns of D ' z h s(ld ' (y^; //, 0, ^) 
and D z h s<ld) (yid, fi,0,v) are indeed equal. Consequently the first row of their adjugates is 
the same: 

e[adj(I - D z h s(ld) (y ld] //, 0, v)) = e]adj(I - D z h S(ld) (y ld ; /i, 0, v)) . (46) 

The combination of ( 1431) . ( l44l) . ( 1451) and = produces 

(l-D z h s(ld \y ld ;v,Q,iy))z ld {fi,0,v) = (i - D z h S(ld) {y ld ; /i, 0, v))z ld {^ 0, v) + 0(3) , (47) 

and multiplication of both sides of (j4"T|) by f|46l) on the left (remembering (Q and = e^^i) 
leads to 

det (/ - D z h s(U) {yi d \ /i, 0, v))uu{^ 0, v) = det (/ - D z h s(U) {y u ; /i, 0, v))ui d (^, 0, z/) + 0(3) . 

We now plug in known expansions (ui d is given by fl3il . det(J — D z h s(ld] ' (yi d ; /i,0, z/)) = 
5 + 0(1) by ([35]) and Lemma El det(I - D z h s ° d) (y ld ; /i, 0, v)) = k fx + k 2 v + 0(2) by ([36]), 
(142|) and Lemma [3]): 

+ fc 2 z/ + 0(2))i/(l + 0(1)) = (5 + 0(l))u u (ji, 0, i/) , 
which, upon rearranging, produces the following useful expression: 

uufr, 0,1/) = hpv + hv 2 + 0(3) . (48) 



Below we use the same approach to derive further expressions for u and u\ d with different 
assumptions on the parameters. For brevity we will not provide the same level of detail. 

b) Suppose i] = and \i = 0. As above, since 77 = 0, the 5-cycle is also an 5-cycle. In 
particular, z_ d is a fixed point of h s( d) . Here we express h s( d) as a Taylor series centered 
at yo (the reason for choosing y will soon become clear) and evaluated at z_ d : 

MO* 0. ") = hS( ~ d) fool 0, 0, u) + D z h s( ' d) (y ; 0, 0, u) (z_ d (0, 0, 1/) - y ) , (49) 

where, unlike in (j4"3"j) . there is no next term in the expansion when /x = the map is affine. 
Now, zq is a fixed point of h s , therefore 

£b(0, 0, 1/) = ^(y ; 0, 0, v) + D z h S (y - 0, 0, v) {z (0, 0, 1/) - y ) . (50) 

But S^~ d ^ = <S® (by ( |T2|) ). and since to = (Lemma EKE])) we may perform the same simpli- 
fication that we did above to combine f l4"9]) and ( 150|) leaving 

(l-P z /i 5C ~%o;0,0,i/))MM,*') = (/-^(yo;0,0,z/))i (0,0,z/) , 
det(/-^/i 5C " d) (yo;0,0,i/))M_ rf (0,0,^ = det (/ - D z h S (y ; 0, 0, i/))uo(0, 0, 1/) , 
(A; 2 z/ + 0(z/ 2 ))(t_ d + 0(z/)) = (5 + O(z/))n (0,0,j/), 

=> £o(0,0,i/) = ^t^z/ + 0(z/ 2 ) . (51) 
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c) Suppose v = 0. As discussed in §|6j when v = the 5-cycle is also an d )-cycle. 
Therefore Zq is a fixed point of h s( d) and so 

zobi, 77, 0) = h s '~ d) (y ; /x, 77, 0) + D z h s( ~ d) {y ; /x, 77, 0) (z {fx, 77, 0) - y ) +0(3). (52) 

Also, Zo is a fixed point of h s , thus 

5 (/i, ?7, 0) = ^(y ; 7/, 0) + D z h S (y ; /x, rj, 0) {zq(ji, tj, 0) - y ) + 0(3) . (53) 

Combining ( 152"]) and fl53|) leads to (since «S^ - ^ = «S°) 

det (/ - Z^/i 5 ' d> (y ; A*, »7, 0))tto(//, 7/, 0) = det (/ - D z h s (y ; /x, 77, 0)) w (//, 77, 0) + 0(3) . 
Take care to note that 

fco = |- det(J - D x f(x d (0; 0) , (54) 
a/x C=o 

is different to fl32|) due to the presence of nonlinear terms in fl6]) (below we will show that in 
fact k = —ko). Consequently 

(k 0f x + k lV + O(2))rj(l + O(l)) = (5 + O(l))u (n,v,0) , 

«o(^r/,0) = ^+^ 2 + 0(3). (55) 
o o 

d) Suppose v — and /x = 0. Here -Z(j+i)d is a fixed point of /x 5<w) , so 

%+i )d (0, 77, 0) = h sild) {y ld ; 0, 7?, 0) + D z h s ™ {y ld ; 0, 77, 0) (i (J+1)d (0, 77, 0) - y ld ) , (56) 
and Zia is a fixed point of h s{ld) , so 

l w (0, 77, 0) = h S(ld) {y ld - 0, r/, 0) + D z h m) {y ld ; 0, 77, 0) (z ld (0, 77, 0) - y ld ) , (57) 

and since 

S (id) = S (id)o we obtain 

det (/ - D z h sm (y ld ; 0, 77, 0))« ({+1)d (0, 77, 0) = det (J - D z h^° d) (y ld ; 0, 77, 0))« w (0, 77, 0) , 
(A; 1 r7 + 0(r7 2 ))(t (m)d + 0(7 ? )) = (5 + 0(77)^(0,77,0), 

u w (0, 77, 0) = ^%^77 + 0(77 2 ). (58) 



We now apply the implicit function theorem to the C K ~ 1 function Ui d (^). By (f48|) and 
( )58l . there exists a unique C K ~ l function <px such that for small \i and z/, w/d(/x, 0i(/x, ^), v) = 
and 

<M/x, 1/) = -j^—liu - T p—u l + 0(3) . (59) 

Recall that ui d = along the /x-axis (a consequence of fl30|) and fl3Tl) ). therefore 0i = 
whenever v — 0. Thus we may rewrite f l59|) as 



= ;/(-— ^/x-—^—^ + 0(2)) . (00) 

Kit(i+i) d 
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Similarly by (|5Tj) and ( )55|) . there exists a unique 1 function 2 such that for small /i and 

77, u (n,r],(f) 2 (ix,r])) = and 

2 (/i, V ) = V (--Jj-fJi ~ t^V + 0(2)) , (61) 
^ k 2 t- d k 2 t- d J 

where the r\ may be factored in the same fashion as for ( |60l) . 

Step 2: Verify (TJj o/ £/ie lemma and refine the formulas for <pi and <p 2 , ( TffOj) and ( fgij) . 
Here we continue to employ the methodology above to compare the S and 5-cycles. 

a) Suppose r\ = (f>i(fi, v). Then = 0, and so the 5-cycle is also an iS-cycle. Thus, in 
particular, zq is a fixed point of h s : 

+ D z h s (y Q \^(j)x{^ u),u)(zo(fM,(pi(li, u),u) - y ) +0(3) . 
Also Zq is a fixed point of h . 

zo{vAi{v,v),v)) = h (y ;ii,(j)i(fM,u),u) 

+ D z h s (y Q ;n,(j)x(n : v), v) (z (//, </>i(/i, ^),^) - y ) +0(3) . 

Then, since S = S°, 

det (/ - D z h s (y ; fi, 0x(/i, u), u)))u (v, ^), i/)) = 
det (/ - D z h S (y ; /i, 0i(/i, v), u)))u (fi, 0i(/i, v), v)) + 0(3) . (62) 

Unlike for similar expressions in Step 1, we have already determined an expansion for each 
of the four components of f[6"2"j) (in particular Uo(fi, (fii(n, v), v) is given by floT]) . ([So"]) and ([6"Uj) 
and Mo(/-t, 0i(/ i ; z/ ) ; ^) is given by ( 13"U|) and fIBTJ]) . also recall Lemma E]): 

(V + fc 2 f + 0(2)) {^f±v + 0(2)) = 

(* + <>{i)) (hs^ - + °( 3 )) + °( 3 ) • 

Equating the second-order coefficients produces 

kiht(i+i)dt-d = -55 . (63) 

b) Suppose r] = <pi(fi,h / ) and /i = 0. In the same manner as above, equating the first 
components of the power series of h s{{ ' 1)d) and h yields 

det(J - D z h sW ~ 1)d) (y ld - 0, 0!(O, u), i/)« ( ,_ 1)d (0, 0i(O, 1/), 1/) = 
det(J - D z h m (y ld ; 0, &(0, */), i/)u w (0, &(0, i/), !/) , 

and since S^" 1 )*) = <S^°, 

(k 2 u + 0(u 2 ))(t^ l)d + 0(u)) = (5 + 0(u))(u + 0(u 2 )) , 

=> ^2 = . (64) 



t{i-i)d 
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c) Suppose v = (j) 2 ((j,,r}) and \i = 0. Here, equating the first components of h s and h s 
produces 

det(J - D z h s {y ; 0, V , 2 (O, rj))u d (0, V , 2 (O, V )) = 
det(J - D z h S (y 0] 0, 77, 2 (O, 7?))u (0, tj, 2 (O, 17)) , 



from which it follows that 

(k lV + 0( V 2 ))(t d + 0( V )) 
=>- ki 



(5 + 0( V ))( V + 0( V 2 )) 



(65) 



We now combine above equations to demonstrate some parts of the lemma. By combining 
(|63|) . (164]) and (1631) we obtain (0). Combining ( )60|) . ( 16^|) and (165]) verifies (Jn|). Combining 
(15T1) . (I6"lj) and (1631) will verify (Jul]) once it is shown that k = —k (see below). 

Step 3: Derive and analyze the one-dimensional center manifold of f s to obtain the 
function A(£) and identify saddle-node bifurcations of S-cycles. 

When £ = (//, 77, i>) = 0, x = is a fixed point of f s (x;t;) and the associated stability 
multipliers are the eigenvalues of D x f s (0; 0) = M s (0). The matrix M S (Q) has an eigenvalue 
1 of algebraic multiplicity one. Let v G M N be the associated eigenvector, i.e. Ms(0)v = v. 
Notice v ^ implies ejv ^ since if not then M s o(0)v = v which contradicts the assumption 
that I — M§(0) is nonsingular (Def. [3]). In what follows we assume ejv = 1. 

We now compute the restriction of f s to the one- dimensional center manifold. Let 



denote the (JV + 3)-dimensional, C , extended map. The Jacobian, 



DF(0;0) 



has a four-dimensional center space, _E C G M^" 1 " 3 , spanned by 



M s (0) 


Ps(0)6(0) 








(66) 



V 




yo 

















l 






















1 
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since, in particular, y = h s (y ; 0) = Ms(0)yo + Ps (0)6(0). 

Since ejv = 1, we may use the center manifold theorem to express the local center 
manifold, W c , of f s , in terms of s = ejx and £. In particular, on W c , 



x 



X{s;Z) = sv + fiy + O{2) 



(67) 
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where X : K x M 3 — > M N is C K 1 . The first component of the restriction of f s to W c is 
given by 

s'(s;0 = e[/ 5 (X( S ;0;0 

= eJ>P 5 (0)&(0) + e[M s (0)(sv + m ) + 0(2) 

= eJ(P s (0)b(0) + M s (0)y o )fi + ejM s (0)vs + 0(2) 

= s + 0(2), 

since e](P s (0)b(0) + M s (0)y ) = ejy = t = and ejM s (0)v = ejv = 1. However we 
require the knowledge of second order terms of s'(s; £), so write 

s'(s;£) = s + cis 2 + C2/is + c 3 r7s + c 4 i/s + C5/i 2 + c 6 /ir7 + C7/iz/ + C8?7 2 + C9?7z/ + cioz/ 2 + 0(3) . (68) 

We now utilize known properties of / to determine an expression for the majority of the 
coefficients, q. 

1) When /i = 0, x = is a fixed point of / 5 , thus s'(0; 0, r], v) = 0, hence eg = eg = cio = 0. 

2) When rj = 0, x = xq is a fixed point of / 5 and since here s = sq = 0, we have 
s'(0; /i, 0, v) = 0, hence 05 = 07 = 0. 

3) Similarly when v — 0, fixed point of / 5 . Here s = s d = u d \i = t d fi + 0(2), 
thus s'(iidfi] fi, T], 0) = iidfM, leading to 

C2 

ci = -— , c 6 = -c 3 t rf . 
td 

Thus we have reduced the center manifold map fl68|) to 

s'(s; = s - ^s 2 + c 2 /is + c 3 77s + Cavs - c 3 t d jj,7] + 0(3) . (69) 
td 

4) Let A(?7, v) be the eigenvalue of D x f s (0; 0, 77, 1/), O x_1 dependent on 77 and z/, for which 
A(0, 0) = 1. This eigenvalue is also the stability multiplier of the fixed point, s = 0, of 
s'(s; 0, 77, u), that is, by (1691 

cV 

^ = -5-(°5 0, 77, 1/) = 1 + c 3 r/ + c 4 z/ + 0(2) . (70) 
as 

Now, det(J — D x f s (0;0,r),v)) is equal to the product of all N eigenvalues (counting 
algebraic multiplicity) of I — D x f s (0; 0, 77, v). Thus 

det(J - D x f(0; 0, 77, 1/)) = (1 - A (77, f))^, , 

where "P(77, z/) denotes the product of the remaining iV— 1 eigenvalues of I—D x f s (0; 0, 77, v). 
V(r), v) is O x_1 and V(0, 0) 7^ since the algebraic multiplicity of the eigenvalue 1 of 
M 5 (0) is one. Let 

k = 7>(0,0), 
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then, using (1701) . 

det{I-D x f{0;0,r),v)) = {-c 3 r) - c*u + 0{2)){k + O(l)) 

= — c 3 k^ — c 4 kz/ + 0(2) . 

Using ( H2l . (101]) and (1051) we arrive at 

5 6 
C3 = — , c 4 



«*(i-i)d 

5) When rj = u = 0, the 5-cycle has two points on the switching manifold and coincides 
with the 5-cycle. Both x (fi, 0, 0) and Xd(fJ>, 0, 0) are fixed points of f s . Let Xi(fi) and 
A 2 (/i) be the respective eigenvalues of D x f s (x (fi, 0, 0); //, 0, 0) and D x f s (xd{(J>, 0, 0); /x, 0, 0), 
O k_1 dependent on jj, with Ai(0) = A2(0) = 1. As above, since so(/x, 0,0) = 0, from 

(EHD 

A x (//) = ^(sqO", 0, 0); fi, 0, 0) = 1 + pa/i + OGu 2 )) , 
and then since Sd(/i, 0, 0) = i^/i + 0(/i 2 ), 

A 2 (//) = — (s d (//,0,0);//,0,0) = l--^(t^ + 0( yU 2 )) + c 2 /i + 0(/. 2 ) 
9s i d 

= l-c 2 /i + 0(/i 2 )) . 

Again, as above, 

det(J-D*/ 5 (z o (/*,0,0);^0,0) = -c 2 Kfi + 0(/i 2 ) , 
and det(I — D x f (x d (/i, 0, 0); 0, 0) = c 2 K/i + 0(/i 2 ) . 

But, recall (1501) . so 

c 2 = , 

and by ([Ml) 

k = -h ■ (71) 

Substitution of (I7T|) into ( 10T|) verifies (Jul]) of the lemma (using also (jMJ) and (1051)). It 
now only remains to demonstrate (jrvj) and (jvj) of the lemma. 

We have shown that the restriction of f s to W c is 

s '( s; £) = s + — s 2 - — /xs r?s us + -fir) + 0(3) . (72) 

Kt d K Kt d K 

We now look for saddle-node bifurcations of ([JJ). Since f£(0;0) = 1 and f^(0;0) ^ 0, 
the implicit function theorem implies that there exists a unique C K ~ 2 function, ip such that 
= 1 for small £ and 

m = jV+£rV+ ^r^v + 0(2) . (73) 
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Then saddle-node bifurcations occur when s'^ (£)>£) = VKO- Let 



A(0 = -^(^(0; 0-^(0)- (74) 

Substitution of f!72|) and ( 175]) into f!74|) yields ( |3"7|) . To complete verification of (jiv]) of the 
lemma we formally show that (1721) has a saddle-node bifurcation at A(£) = whenever fi > 
by verifying that all nondegeneracy conditions of the saddle-node bifurcation theorem [33] 
are indeed satisfied: 

ds' 

1) by construction, — — (/0(Oj£) = 1 when A(£) = 0, 

as 

ds^ St> 

2) — — = — u + 0(2) when it > verifying transversality (where v is given in ( 1401) ). 
ov 2k 

Step 4: Compute the intersections of A = roi/i 77 = and z/ = </> 2 (/i, 77) to obtain the 
functions £i and 

To determine where A(£) = intersects rj = it is natural to look at A(/z,0, v) = 0, but 
this is insufficient for a derivation of Ci(aO because we may not apply the implicit function 
theorem to A(£) since it contains no linear terms. Instead we use the fact that when rj = 0, 
the (S-cycle is also an 5-cycle. We have A = when, in addition, this periodic solution has 
an associated multiplier of 1 as an S-cycle. This occurs when (using ( 136|) . (j4"2"|) and (164|) ) 



det(J - D x f*(x (^ 0, i/); //, 0, 1/) = fc /x + v + 0(2) = , 

t(i-i)d 

Application of the implicit function theorem to the previous equation produces 



Moreover, 

AO*, 0, u) = -^—{v - d(/i)) 2 + o(\v- Ci(^)| 2 ) , 

and so for \i > 0, A > on the z/-axis. 

The curve (/i, C2O"), fait^i C2 (/•*)) along which A(£) = intersects the surface 2 (/i,n), is 
easily computed in a similar fashion. When \x > 0, A > along z/ = 2 (yU, 77) an d so A < 
only when 77 < and v > ^(a*, 77) and stated in the final part of the lemma. □ 

Proof of Thm. [73 We begin by determining the region of admissibility of the S-cycle, {x{\. 
From (|28|) . = /i(uj+0(l)), thus since by assumption y and t/ w are the only points of {i/i} 
that lie on the switching manifold for small £ with \i > 0, here lies on the same side of 
the switching manifold as yi for each i 7^ 0, Zd. The n-cycle, {yi}, is admissible by assumption 
(for fj, > 0), thus {£j(£)} is admissible exactly when s , s/d > (since So = Sw = R). By 
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(|30|) and (I3TT) s , sjy > when 77, 1/ > 0, therefore the iS-cycle is admissible in ^ as stated 
in the theorem. 

Similarly, for /i > the iS-cycle is admissible exactly when sq, < (since Sq = S^ = L). 
By Lemma EtHID, s w (f) = when 77 = fa(jJi,v). By ([58]), ([65]) and Lemma 1®, ^k(f) = 

— + 0(£) which is positive for small £, hence Sid < for 77 < 0i(/j, z/) when /j > 0. 
Similarly by Lemma [9 ljm]) . so(£) = when v = 02 (/^, v) and by ([5T]) . (164)) and LemmaHl, 
— — — I" ^(O which is positive for small £, hence Sq < for z/ < 2 (/^ ) ''7) when 

/ii > 0. Therefore the 5-cycle is admissible in ty 2 - 

It remains to verify admissibility of 5-cycles. For the theorem in (TO] this was straight- 
forward since, if it existed, the iS-cycle was unique. The situation here is more complicated 
because there may be two coexisting, admissible 5-cycles. In the proof of Lemma [9] we 
determined the restriction of f s to the center manifold through (s; £) = (0; 0), (172]) . When 
/i > and A(£) > 0, locally the map f[72]) has two distinct fixed points, say, so,i and so,2- 
We will denote the corresponding iS-cycles of ([6]) by {a^i} and {2^,2} and assume so,i > so,2- 
For small £ within {£ | \i > 0, A(£) > 0}, so,i and so,2 are C K ~ l functions of £. 

On the surface A(£) = the two solutions coincide (see ( 1731) ): 

s ,i(0 = so, 2 (£) = when A(£) = 0. 

At the intersection of A(£) = and r\ = 0, namely (/i, 0, Ci(aO) ( see Lemma [9](|vj)), so,i = 
s , 2 = 0. Thus by (1751) . on A(£) = 0, s 0) i = £0,2 < when z/ < Ci(^) and s ,i = s ,2 > when 
z/ > CiO")- 

Now, for fi > 0, So,i and so )2 can only be zero if 77 = because if s = is a fixed point of 
(1721) then the corresponding 5-cycle would also be an 5-cycle which for /j > must be {ii}, 
so then so(0 = and by (1501) we would necessarily have 77 = 0. Consequently one of so,i 
and so, 2 is zero when rj = and so,i and $0,2 are both nonzero when 77 7^ 0. Since we assume 
so,i > so t 2 for A(£) 7^ 0, when /n > and 77 = we must have so,i = when v < Ci(yu) and 
s 0) 2 = when v > Ci(aO- Consequently s ,i < in U ^3 and so,2 < in ^1 U \&2 U ^3. 

Before we are able to perform a similar analysis of for i 7^ 0, we find it necessary to 
first derive an expression for sy in terms of ti and t^+i- Recall that the center manifold, 
W c , is given by (1671) where Ms(0)t> = v and e]"f = 1. When £ = 0, t/q and y<j are both fixed 
points of h s , thus (I - M s {0))y = P 5 (0)6(0) = (I - M 5 (0))y d and so 

(I-M s (0))(yo-Vd) = 0, 

and since 7/0 7^ 2/d (Lemma [7]) and the eigenvalue 1 of the matrix Ms(0) has algebraic 
multiplicity one, t/ — y d is a scalar multiple of tj. Due to the specified vector scaling we have 

v = ^{y d -yo)- (75) 

Combining ( 157)) and (175)) yields 

= (a* - 2/0 + + 0(2) . (76) 
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This may be generalized to an expression for x it j(£) using x i+ ij = fib + A Si Xij + 0(2) and 
y i+ i = //6(0) + A s .(0)yi, from which we deduce 

SiAO = (m - ^) U + + 0(2) , (77) 

and hence 

- s^(0 = -^fc - W(so,i(0 - s 0l2 (0) + 0(2) . 

Therefore for small £ > with > and A(£) > 0, s^i > s i)2 if tj > td+i and s^i < Sj )2 if 
U < td+i (because we assumed so,i > so, 2 ). 

Above we showed that when rj = v = and \i > 0, so,i = and so,2 < 0. Thus here 
sid,i — and s^,2 > and by Lemma [91 along 77 = 0i(/i, f), s^i = and s^ )2 > 0. From this 
is easily follows that s^,i > in \l/ 2 U \l/ 3 and s^, 2 > in ^ U \l/ 2 U \&3. Similarly when v = 
and /i > 0, S(z_i)d,i < and S(|_x)d,2 = and consequently sq^i^i < in U \l/ 2 U ^3 and 
s (i-i)d,2 < in \&x U ^3. By analogous arguments, s_ d l > in ^ U \I/ 2 U ^3 and s_d,2 > in 
^1 U \J/ 3 . By (1771) . for i ^ 0,(1 — l)d, Id, —d, s^i and s i)2 have the desired sign for admissibility 
for small £ with fi > 0. The above statements show that {2^1} is admissible in \l/ 2 U ^3 and 
{xj j2 } is admissible in ^ U ^3 which completes the proof. □ 
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